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Abstract 

< 

! An overview of the cosmological models with expansion, shear and rotation 

, is presented. Problems of the rotating models are discussed, their general 

, kinematic properties and dynamical realizations are described. A particular 

I attention is paid to the possible observational effects which could give a defi- 

' nite answer to the question: is there evidence for the rotation of the universe? 

O ■ 1. INTRODUCTION 

W 



In his lecture "Spin in the Universe" |[432|| , E.T. Whittaker has drawn attention to the 



^ ' following questions: "Rotation is a universal phenomenon; the earth and all the other mem- 



bers of the solar system rotate on their axes, the satellites revolve round the planets, the 
planets revolve round the Sun, and the Sun himself is a member of the galaxy or Milky Way 
system which revolves in a very remarkable way. How did all these rotary motions come into 
being? What secures their permanence or brings about their modifications? And what part 
do they play in the system of the world?" 

The standard Friedman-Lemaitre cosmological models are based on very idealized as- 
sumptions about the spacetime structure and matter source. The Friedman-Robertson- 
Walker metrics describe an exactly isotropic and homogeneous world filled with an ideal 
fluid (radiation and dust are the well known particular cases). Although these assumptions 
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are qualitatively supported by observations, it seems doubtful that the geometry of the uni- 
verse and its physical contents are so finely tuned starting with the very Big Bang till the 
modern stage of the evolution. A reasonable step is to look for more general and realis- 
tic models with the assumptions of isotropy, homogeneity, and ideal fluid matter contents, 
relaxed. Here we will not discuss inhomogeneous cosmologies which are reviewed, e.g., in 



m 



Already in the early years of the general relativity theory, attempts were made to con- 
struct (exact or approximate) solutions of the gravitational field equations with rotating 



matter sources (see, e.g. ||225| , pT|jT^ , |23 1| , |407|| ) . Lanczos |p25|| appears to be the first who con 



sidered the universe as the largest possible rotating physical system in a model of a rigidly 
rotating dust cylinder of an infinite radius. Dust density in this solution (which was later 



rederived by van Stockum [|407|| ) diverges at radial infinity, which was a serious problem. 



In 1946, Gamow ||112|| stressed a possible significance of the cosmic rotation for explaining 



the rotation of galaxies in the theories of galaxy formation, see also his discussion in ||113 



Independently, Weizsacker [ [43(J| ] initiated the study of the primordial turbulence and of 



its role in the formation of the universe's structure. Soon after this, K. Godel | p.l6| | had 
proposed a stationary cosmological model with rotation in the form (coordinates {t, x, y, z}, 
a =const): 

ds^ = (^dt^ - 2e''dtdy + ^e^'^dy^ - dx^ - dz^^ . (1.1) 

In the Godel's model, the matter is described as dust with the energy density e, and the 
cosmological constant A is nontrivial and negative (i.e. its sign is opposite to that intro- 



duced by Einstein). The angular velocity uj of the cosmic rotation in (|1 . 1| ) is given by 
cu^ = 2^ = AnGe = —A, with G as Newton's gravitational constant. For many years 
this model became a theoretical "laboratory" for the study of rotating cosmologies, see 
e.g. [|6ip9|,p^,[r9|, [l5^ , |l5^jm|j336H338^ . The stationary solution ([l^) is distinguished 



among other spatially homogeneous models: the so-called Godel theorem (complete proof is 
given in ||300|| ) states that the only spatially homogeneous models with zero expansion and 



shear are Einstein's static universe and Godel model (|1 . 1| ) . 

Subsequently, Godel ||117|| himself outlined a more physical expanding generalizations of 



( p..lD , although without giving explicit solutions. Later a considerable number of exact and 
approximate models with rotation and with or without expansion were developed. At the 
end of this paper one can find a list of references which may be considered as a replacement 
of a detailed review of the previous work on rotating cosmology. 



Here we would like to mention only briefly some earlier papers. Maitra [p42|] gave a 
generalization of van Stockum's solution which is cylindrically symmetric inhomogeneous 
stationary dust-filled world with rotation and shear. Maitra formulated a simple criterion 



of absence of the closed timelike curves. Wright |[437|| presented another inhomogeneous 



cylindrically symmetric solution for dust plus cosmological term. 

Dust and cosmological constant represented the matter source in the series of papers by 



Schiicking and Ozsvath [P95| - p04| , where the spatially homogeneous models with expansion, 
rotation and shear were considered. The closed Bianchi type IX solution ( ||299|| with sta- 
tionary limit | |295| , |296|| ) has drawn a special attention. Most general anisotropic Bianchi type 
IX model was studied by Matzner | p^61| , p^62| ] who took ideal fluid plus collisionless massless 
particles (photons and neutrinos) as matter source. The approximate solution was given. 
Later, Rebougas and de Lima [p59|] described the exact analytic Bianchi type IX solution 
with expansion, shear and rotation for a source represented by dust with heat flow and cos- 



mological constant. Analogously, Sviestins [[412|| presented an exact Bianchi type IX solution 



for ideal fluid plus heat flow. In an interesting recent paper, Gr0n [|127|| gave a shear-free 
Bianchi type IX solution in absence of matter (but with cosmological constant) which can be 
considered as a rotating generalization of de Sitter spacetime. In fact, the latter is recovered 
asymptotically after the quick decay of the vorticity. 

The wide spectrum of spatially homogeneous models includes approximate and exact 
solutions for all Bianchi types and different matter contents (from simple dust to sometimes 
quite exotic energy-momentum currents). To mention but a few, Ruzmaikina-Ruzmaikin 
|375|| , Demiahski-Grishchuk [^, and Batakis-Cohen presented approximate solutions 



of Bianchi type V and VII with expansion, shear and vorticity for the ideal fluid in non- 



comoving coordinates. Bradley and Sviestins |^ added heat flow to the ideal fluid source 
in their analysis of a Bianchi type VIII model. In a series of papers, Fennelly ||104| - |108( 
has studied Bianchi models of different types for homogeneous universes filled with perfect 
charged (conducting) fluid and magnetic field. A neutral ideal fluid plus electromagnetic 
field and null radiation represented the matter source in a shear-free expanding and rotating 
solution of Patel-Vaidya [[42 2|] which can be interpreted as a standard Friedman-Robertson- 
Walker cosmology with rotating matter [|130|| . More sophisticated matter source may include 
viscous fluid and scalar fields, see, e.g., recent papers |p46H256|| . More exact solutions, which 
we will not discuss in detail, M,p6iP, pO| , p83| , Plj356|j357| , |370| , [412| ,^ were reported in the 



literature. Most of these models were overviewed by Krasihski ||198H202|| , more recent review 
of generalized cosmologies with shear and rotation see in p06| - |208| . 



It is worthwhile to notice that along with the constant and deep interest in the rotating 
cosmologies, historical development has revealed several difficulties which were considered 
by the majority of relativists as the arguments against the models with nontrivial cosmic 
rotation. One of the problems was the the stationarity of Godel's model ( |1.1| ). Apparent 
expansion of the universe is usually related to the fact of the red shift in the spectra of distant 
galaxies, and thus all the standard cosmological models are necessarily non-stationary. It 
was discovered, however, that it is impossible to combine pure rotation and expansion in a 
solution of the general relativity field equations for a simple physical matter source, such as 
a perfect fiuid [p5H 97[|. There are two ways to overcome that difficulty: one should either 
take a more general energy-momentum or to add cosmic shear, the explicit examples were 
already mentioned above. The possibility of combining cosmic rotation with expansion was 
the first successful step towards a realistic cosmology. 

Another problem which was discovered already by Godel himself, was related to causality: 
the spacetime ( |1.1| ) admits closed timelike curves. This was immediately recognized as 
an unphysical property because the causality may be violated in such a spacetime, see 
relevant discussion in [|118| , |143|,[l34l , |T35| , p27|J336lj406| . Considerable efforts were thus focused 



on deriving completely causal rotating cosmologies. In his last work devoted to rotation, 
Godel without proof mentions the possibility of positive solution of the causality problem 



117|| . First explicit solutions were reported later ||242| , |297|| . Maitra ||242|| formulated a simple 
criterion for the (non)existence of closed timelike curves in rotating metrics which we will 
discuss below. Eventually, it became clear that the breakdown of causality is not inevitable 
consequence of the cosmic vorticity, and many of the physically interesting rotating models 
were demonstrated to be completely causal. That was the second important step in the 
development of the subject. 

The discovery of the microwave background radiation (MBR) has revealed the remarkable 
fact that its temperature distribution is isotropic to a very high degree. This fact was for a 
long time considered as a serious argument for isotropic cosmological models and was used for 
obtaining estimates on the possible anisotropics which could take place on the early stages of 
the universe's evolution. In particular, homogeneous anisotropic rotating cosmologies were 



analyzed in |p.42| , pl| , p3| , p^ , |18CI|J435|| , and very strong upper limits on the value of the cosmic 
rotation were reported. Unfortunately, usually in these studies the effects of vorticity were 
not separated properly form the cosmic shear .[| For the first time, in [ [^87| , ^88| ], the effects 
of pure cosmic rotation were carefully separated from shear effects. The most significant 
result is the demonstration that in a large class of spatially homogeneous rotating models 
the pure rotation leads neither to an anisotropy of the MBR temperature nor to a parallax 
effects, although the geometry remains completely causal. A typical representative of that 



class is the so called Godel type metric with rotation and expansion ( [4.1|) which is considered 
below. Most recent analyzes, based also on the CORE data on the anisotropy of MBR, see 
|6T| , |237| , p!80| , again estimated not a pure rotation but rather a shear. 



The last but not least problem is the lack of direct observational evidence for the cosmic 



^Even now in the literature one can often encounter statements like: "If the universe had vorticity, 
the microwave background would be anisotropic..." [[403||. 
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rotation. Attempting at its experimental discovery one should study possible systematic 
irregularities in angular (ideally, over the whole celestial sphere) distributions of visible 
physical properties of sources located at cosmological distances. Unfortunately, although 
a lot of data is already potentially accumulated in various astrophysical catalogues, no 
convincing analysis was made in a search of the global cosmic rotation. Partially this can 
be explained by the insufficient theoretical study of observational cosmology with rotation. 
To our knowledge, till the recent time there were only few theoretical predictions concerning 
the possible manifestations of the cosmic rotation. Mainly these were the above mentioned 
estimates of MBR anisotropics [|142| , p| , p^ , |180|| (also of the X-ray background anisotropics 
435|| ) and the number counts analyses ||225| , |1 17| , [431| , P65| , |103|] . Besides that, it was pointed 



out in pi7| , p73|| that, in general, parallax effects can serve as a critical observational test for 
the cosmic rotation. 

Very few purely empirical analyses (without constructing general relativistic models) of 
the angular distributions of astrophysical data are available which interpreted the observed 
systematic irregularities as the possible effects of rotation. Of the early studies, let us 
mention the reports of Mandzhos and Tel'nyuk-Adamchuk | p^43| - |245| on the preferential 
orientation of clusters of galaxies, Valdes et al | |424| ] on the search of coherent ellipticity of 
galaxy images, and Andreasyan ||l6|,|l7l on flattening of distant galaxies. 



Birch P6| , |47| , |82| (Jodrell Bank Observatory) reported on the apparent anisotropy of the 
distribution of the observed angle between polarization vector and position of the major axis 
of radio sources, and related this global effect (using heuristic arguments) to the nontrivial 
cosmic rotation. Subsequently, these observations were disputed from a statistical theory 
viewpoint P35| , ^^ , but the analysis using indirectional statistics ||173|| confirmed Birch's 
results. It seems worthwhile to mention that the independently data from the Byurakan 
radio observatory indicated the same anisotropy along, approximately, the same di- 

rection. These observations have stimulated the interest in rotating models; nevertheless, 
no further observations were reported until the most recently, when the new data appeared 



2 79| , |2 79| , p47| , p48|| claiming the same effect. Note however, that again the statistical signifi- 
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The theoretical analysis of the observational effects in the class of viable rotating cos- 
mologies ||287| , |187| , |288| , |194|] has shown that the cosmic rotation produces a typical dipole 
anisotropy effect on the polarization of electromagnetic waves. Considering the case of the 
distant radio sources, one finds the dipole distribution |p.87| , p.94|| r] = uqt cos 9, see ( [7.32[ ) 
below, for the angle rj between the vector of polarization and an observed direction of the 
major axis of the image of a radio source. Here uq is the present value of the cosmic ro- 
tation, r is the (apparent area) distance to a source, and 6 is the spherical angle between 
the rotation axis and the direction to the source, Z is the red shift. The data of Birch and 



the new data of Nodland and Ralston [P79| , P79| , P47|J348|| may be extremely important for 
improving the direct estimates of vorticity. 

The study of observational effects in rotating and expanding cosmological models is 
essentially based on the kinematical properties of spacetime. The next step is to investigate 
the dynamical aspects of cosmological models with rotation, in other words, to construct 
viable models as solutions of the gravitational field equations. The ultimate goal is to 
achieve a self-consistent description of the physical matter sources and to establish a complete 
scenario of a realistic evolution of the universe, including the most early stages of the big 
bang. Already Lemaitre [p30|| and Whittaker [|432|| raised questions about the nature of 
a rotating "primeval atom" from which the universe is supposed to emerge. A related 
problem is the developments of matter irregularities which is relevant to the galaxy and large- 
scale structure formation. Gamov [|112| , |113| , p!^ and Weizsacker |[430|| stressed the importance 
of the cosmic vorticity in the early universe. Silk ||39iy392|| was the first to demonstrate 
instability of rotating models to axial density perturbations and stability to perturbations 
in the plane of rotation. The idea that primordial vortical motions of cosmological matter 
gave rise to the formation of large-scale structure underlied the related study of the so-called 
whirl theory (cosmological turbulence) developed some time ago by Ozernoy and Chernin 
|7|-^^|Tg,gT|-g53],g, see also p75|PB| ,|22|P|,|g^, pB| , pg| , |5g3i . In the recent paper Li 
232|| has demonstrated that the global cosmic rotation can provide a natural origin of the 



rotation of galaxies, and the empirical relation J ~ M^/^ between the angular momentum 
J and the mass M of galaxies can be derived from the cosmological vorticity. [See Brosche 
57|-|60| for empirical data, and Tassie |[413|| for the string-theoretical discussion]. 



Of the recent theoretical developments, it is worthwhile to mention the work [123 
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which Grishchuk has demonstrated that the rotational cosmological perturbation can be 
naturally generated in the early universe through the quantum-mechanical "pumping" mech- 
anism similar to the primordial gravitational waves. These perturbations were considered as 
the source of the large-angular-scale anisotropy of MBR, and comparison with observations 
was made. However, as was recognized by Grishchuk, the pumping mechanism only works 
when the primeval cosmological medium can sustain torque oscillations. No attempt was 
made to study the possible physical nature of the relevant cosmological matter. Exploit- 
ing the close relation between rotation and spin, one can try to attack this problem using 
the models of continua with microstructure. Self-consistent variational general relativistic 



theory of spin fluids was constructed in ||289| - p9l[] . The investigation of the development of 
small perturbations for the models with spinning cosmological matter ||29 1|| yielded a gener- 



alization and refinement of the result by Silk ||392|| about instability to perturbations along 
the spin. 

Most of the dynamical realizations of rotating cosmological models were obtained in 
the framework of general relativity theory as the exact or approximate solutions of Ein- 
stein's gravitational field equations in four dimensions. However, we can mention also 
the study of higher- derivative extensions of general relativity [|I|-^ |157[| , solutions in 5 di- 
mensions [|218|J365|[ , in 3 dimensions [[L33| , [433|J434|[ , and in string-motivated effective cosmol- 



ogy [p3| , pT| . On the other hand, rotation, spin and torsion are closely interrelated in the 
Poincare gauge theory of gravity |[L48| , p.49| , p41|J168|J266|| , and hence it is quite natural to 
study cosmologies with rotation within the gauge gravity framework. Exact rotating so- 
lutions in Einstein-Cartan theory of gravity with spin and torsion were described, e.g., in 
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,15, 
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,289|,|308|,P 


09|,p29|J40q,401 



401 [. General preliminary analysis of the separate 
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stages of the universe's evolution was made in our works ||164| - |166(] , while in ||189| , p,67|| the 
complete cosmological scenarios are considered. 

Summarizing the present state of the problem, it is noteworthy once again to quote 



Whittaker ||432|| : "It cannot be said, however, that any of the mathematical- physical theories 
that have been put forward to explain spin (rotation) in the universe has yet won complete 
and universal acceptance; but progress has been so rapid in recent years that it is reasonable 
to hope for a not long-delayed solution of this fundamental problem of cosmology". 

2. STATIONARY COSMOLOGICAL MODEL OF GODEL TYPE 

The Godel metric (|1 . 1| ) represents a particular case of a wider family of stationary cos- 
mological models described by the interval 

ds^ = dt^ - 2^e'^''dtdy - {dx" + ke^'^'^dy^ + dz^), (2.1) 

where m, a, k are constant parameters. Clearly, cr > 0, and for definiteness, we choose m > 0. 
The metric (^TTP is usually called the model with rotation of the Godel type. Coordinate z 
gives the direction of the global rotation, the magnitude of which is constant: 



/i m a , , 

.= y'-.„„.- = -y'^. (2.2) 

As we see, vanishing of m or/and a yields zero vorticity. 

Let us describe the geometry of a Riemannian spacetime in detail. The determinant 
of the metric is equal (detg) = —{k + a)e^'^^ , and hence we assume that {k + a) > in 
order to have the Lorentzian signature. 

Technically, all calculations are simplified greatly when exterior calculus is used. It is 
convenient to choose at any point of the spacetime (|2.1|) a local orthonormal (Lorentz) tetrad 
/i" so that, as usual, g^^ = h'^h'^riab with riab = diag(+l, —1, —1, —1) the standard Minkowski 
metric. This choice is not unique, and we will use the gauge in which 



/iO = 1, hl = -v^e'"^ h\ = hl = 1, hj = e'^'^VkTa. (2.3) 



Hereafter a caret denotes tetrad indices; Latin alphabet is used for the local Lorentz frames, 
a, b, ... = 0, 1, 2, 3. The local Lorentz coframe one-form is defined, correspondingly, by 

r := h^^dx^. (2.4) 

The curvature two-form for ( PTT| ) is computed straightforwardly: 



Jiab ^^2^a ^^b ^___^^a^b_^b^a\^l ^^2 a, 6 = 0, 1, 2. (2.5) 



k + a 

It is worthwhile to note that the spacetime ( p.l|) has a formal structure of a product of the 
curved three-dimensional manifold with the curvature ( p.5|) times fiat one-dimensional space 
(z-axis), which is reflected in the range of indices in ( |2.5|) . 

The only nontrivial components (with respect to the local Lorentz frame (|2.3| )) of the 
Weyl tensor follows: 

^Oioi ~ ^0202 ~ ~ *-^2323 ~ ~ ^3131 ~ 2 ^1212 ~ ~ 2 ^0303 ~ ( ^ _|_ ^ I • (2-6) 



One thus verifies that the model ( p.l| ) belongs to type D according to Petrov's classification. 

The first dynamical realization of the model ( p.l| ) was obtained by Godel ||116|| . In general 
relativity theory, metric satisfies Einstein's field equations: 

Gab '■= Rab — ^dabR = ^ Qab + I^Tab. (2.7) 



Godel ||116|| considered the simplest matter source: ideal dust with the energy-momentum 
tensor Tab = pUaUb. Substituting the four-velocity of the co-moving matter = 5?, and 
using the components of the Einstein tensor, 

Goo = (l + 4^) , Gii = G22 = a;2, G^^ = 00' U + a'^Y (2.8) 



one finds the parameters of the famous Godel solution [compare with ( |1.1| )]: 



k 1 

cj^ = - A, Kp = 2uj^, - = - -. (2.9) 

a 2 

As we see, both the cosmological constant A and the parameter k must be negative. 
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Without confining ourselves to a particular dynamical realization, we should study the 
whole class of stationary models ( |2.1| ). There is a wide symmetry group for that metric. 
Three evident isometries are generated by the Killing vector fields 

e(o) = |, ^(1) = ^, ^(2) = ^. (2.10) 

These are all mutually commuting. However, one can find two more Killing vectors. Most 
easily this can be done with the help of the suitable coordinate transformations. 
Without touching z coordinate, let us transform the rest three 

(t, x, y) — > (r, r, Tp) 

as follows: 



e"^^' = e'^ cosh (mr), (2.11) 
Qh(mr 



ye-='-^^^ (2.12) 



tan 



- sinh((mr). (2-13) 
Straightforward calculation gives a new form of the metric (|2.1|) : 



.2 _ ,^2 sinh(mr) + a] , 2 jV 



ds = dr — dr -— : du? — dz H J- smh.(mr) dr d(p (2.14) 

m^{k + a) m]lk + a ^ ^ ^ v y 

= (dT+-, sinh(mr) #V - f dr' + 1""^^ + dzA . (2.15) 
\ k + a J \ j 

We thus immediately find the fourth Killing vector, dip. In the original coordinates it reads: 

^3 = — = -—-y— (216) 
dip m dx dy 

Alternatively, let us consider another transformation 

(t, X, y) — > (r, r, (f) 

given by the formulae: 
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ye 



cosh (mr) + cosy? sinh (mr) 
snap sinh (mr) 



tan 



m k + a 



it 



my/k + a 
sin (f 



cosy? + coth(f:) 



The metric ( |2.1|) in the new coordinates reads: 

4 sinh^ (f^) [k cosh^ (f:) + a] 



dr' - dr' 



im?{k + 0") 



dip^ — dz^ 



(2.17) 
(2.18) 

(2.19) 



(2.20) 



+ 



my k + a 



sinh 



mr 
. 2 , 



dr d<f 



dr , 



a 



my k + a 



sinh^ [sf) dip 



, 2 sinh (mr) ,2 1 1 
dr^ H — -dif^ + 



(2.21) 



Evidently is also the Killing vector field. In the original coordinates it looks somewhat 

complicate and can be represented as a linear combination of three vectors: 
d 1 



a 



1 



d({) myk + a'"'"' ' 2m^k + a 
where we denoted the fifth Killing vector: 



^(1) - V A; + 0-^(4), 



mx Q d \ 



-Imx 



m(k + o) 



my 



d_ 
dy' 



(2.22) 



(2.23) 



m{k + cr) dt dx 2 
Three of the five Killing vectors (|2.10| ), (|2.16|) , (|2.23|) satisfy commutation relations 

[^(3),'^(i)] = ^(1), [^(3), ^(4)] = -^(4), [^(i),'^(4)] = m^is)- (2.24) 

Thus ^(1), ^(3), 'C(4) form the closed subalgebra isomorphic to so(2, 1). Two Killing vectors 
^(o),^(2) span the center of the isometry algebra, commuting with all its elements. 

In order to get some insight into possible causal problems, it is helpful to compute squares 
of the Killing vectors. Straightforwardly we obtain: 

(e(o)-e(o)) = i, (2.25) 
(e(l)■e(l)) = -A;e2"^^ (2.26) 

(e(2)-e(2)) = -l, (2.27) 



(^(3) -^(3)) = -[ky e 



,2 „2mx 



m^ 



(e(4)-e(4)) = -7le-mr+ . 

4 \ m[k + a) 



(2.28) 
(2.29) 
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The vector field ^(o) is always timelike, its integral curves coincide with the lines of coordinate 
time. The cases k > 0, k < and k = are essentially different, and it is reasonable to 
consider them separately. 

For positive values of k the four Killing vector fields 'C(i)) C(4) are strictly spacelike; 
the first three of them provide spatial homogeneity of the t = const hypersurfaces. We will 
demonstrate the complete causality of the models with k > later. In particular, the closed 
timelike curves are absent in these geometries. 

For negative k, the vectors ^(i) and ^(4) become timelike, whereas ^(3) can be both spacelike 
and timelike, depending on the values of spatial coordinates x, y. One can immediately 
verify the existence of timelike closed curves, which was first noticed already by Godel 
|116|| . Indeed, consider the closed curve {t{ip),x{ip),y{ip)}, with < f < 2tt, given by the 



coordinate transformation ( p.l7| )-( p.l9| ) for r = const, r = const, z = const. From ( p.20|) we 



find that for large enough r, when cosh^(^) > |f |, such a curve is timelike. By construction, 
it is closed, and this may lead to the violation of causality. Although in ||7^ it was shown 
that such curves are not geodesies, one can in principle imagine a motion of a physical 
particle along these curves under the action of certain forces. It is worthwhile to note that 
the form of the metric is similar in two different coordinate systems, compare ( |2.15| ) and 



( p.21| ). However, for A; < the coordinate Tp curve (with r = const, r = const, z = const) is 
also timelike for cosh^ (mr) > |f |, but it is not closed. 

The special case k = was studied in detail by Rebougas and Tiomno |P6(]|| , see also 



||364|| and more recently ||368|| . The Riemannian curvature is simplified greatly for k = 0. In 



particular, the Weyl tensor ( |2.6| ) vanishes, and thus four- dimensional metric is conformally 
fiat. At the same time, ( p.5| ) reduces to R"^^ = uJ^ i?" A#, with a, 6 = 0, 1, 2. The last relation 
describes a three-dimensional manifold of constant curvature, i.e. (anti) de Sitter spacetime. 
This space has much wider isometry group than /c 7^ case. In addition to the 5 Killing 
vectors (|2iq) , p:T6| ), (j^), we easily find 2 more: 



Q e '"^ 

^(5) = sin(mt)— - cos(mt)— + sin(mt) — , (2.30) 
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d d e '"■^ d 

^6) = cos(mt)— + sin(mt)— H ^ cos(mt) — . (2.31) 

ot ox ay 

Together with these vector fields form the second closed subalgebra isomorphic to so(2, 1): 

[^(0),^(5)] ="^^(6), [^(0),^(6)] = -"^^(5), [^(5), ^(6)] = -"^^(0)- (2-32) 



The six Killing vectors ^(a), A = 1, 3, 4, 0, 5, 6 with the commutation relations (|2.24|) , (p.32|) 



describe the complete so(2, 2) algebra of isometrics (algebra of the conformal group of a 
three-dimensional (anti) de Sitter spacetime) [|36^| . It is worthwhile to note that both new 
Killing vectors are spacelike: 

(e(5)-e(5)) = (e(6)-e(6)) = -l, (2.33) 



Starting from the metric in the form ( p.21| ) with /c = inserted, we can make a simple 
coordinate transformation 

(r, r, LP) — > (r, r, Q), 

redefining the angular coordinate, 

Tfl 

^ = -{Q-r\ (2.34) 
which brings the metric to an explicitly (anti) de Sitter form: 

ds^ = cosh^ dr^ - dr^ - sinh^ (sr) dO^ - dz^ . (2.35) 
Note that u; = f for A; = 0. 

3. CLASS OF SHEAR-FREE COSMOLOGICAL MODELS WITH ROTATION 

AND EXPANSION 



In [ p87|J187|| we considered a wide class of viable cosmological models with expansion 
and rotation. Let us describe it here briefiy. Denoting x° = t as the cosmological time and 
X ,2 = 1,2,3 as three spatial coordinates, we write the space-time interval in the form 
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ds^ = dt^ — 2a riidx^dt — o? •jij dx^dx^ , (3-1) 

where a = a{t) is the scale factor, and 

ni = Vaef\ jij = (3abefef\ (3.2) 

Here (a, 6 = 1,2,3) Vail3ab are constant coefficients, while 

e^''^ = ef\x)dx' (3.3) 

are the invariant 1-forms with respect to the action of a three-parameter group of mo- 
tion which is admitted by the space-time (|3.1|). We assume that this group acts simply- 
transitively on the spatial {t = const) hypersurfaces. It is well known that there exist 9 
types of such manifolds, classified according to the Killing vectors C,{a) and their commuta- 
tors [^(a)5^(6)] = C^ab^ic)- Invariant forms ( p.3| ) solve the Lie equations C^^^^^e^""^ = for each 
Bianchi type, so that models ( |3.1D are spatially homogeneous. 

We consider {t, x*) as comoving coordinates for the cosmological matter, that is the four- 
vector of average velocity of material fluid is equal = 6q. It is worthwhile to note that u'^ 
is not orthogonal to hypersurfaces of homogeneity. Such a models are usually called "tilted" . 
The kinematical characteristics of (|3.1|) follows: volume expansion is 

^ = 3-, (3.4) 
a 

nontrivial components of vorticity tensor are 

u^ij = -^C'^ijUk, i,j = 1,2,3, (3.5) 
and shear tensor is trivial, 

a^, = 0. (3.6) 

Hereafter the dot (') denotes derivative with respect to the cosmological time coordinate t. 
Tensor C^ij = e^a)i^i^\'^^ ~ '^j^i'^) is the anholonomity object for the triad (^.31); for I- VII 
Bianchi types values of its components numerically coincide with the corresponding structure 
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constants C°6c- The list of explicit expressions for ^(a),e^°'KC'^bc,C''ij for any Bianchi type 
is given in ||287|| . 



For completeness, let us mention that cosmological models ( |3.1|) have nontrivial acceler- 
ation 

u'^VuU^ = (0, anj ). (3.7) 
We choose the constant matrix I3ab in (|3.2| ) to be positive definite. This important 



condition generalizes results of Maitra ||242|| , and ensures the absence of closed time-like 
curves. 

One can immediately see that space-times ( |3.1| ) admit, besides tree Killing vector fields 
^(a), a nontrivial conformal Killing vector 

Cconi = adt. (3.8) 



All models in the class ( |3.1| ) have a number of common remarkable properties ||287| , |187|| : 
Causality: Space-time manifolds ( |3.1| ) are completely causal. 



Indeed, following p42|| , let us suppose the spacetime (|3rT|) contains a closed curve x'^{s) 



where s is some parameter, < s < 1, which is timelike. The latter means that at 
any point on the curve the length of the velocity vector is strictly positive, 

dx'^ dx" 

Since the curve is closed, x^^{0) = x'^(l), there exists a value Sq at which ^ = 0. At 
this point one has 

dx^ dx" 



ds ds 



= -<.= /3..e.'°'ef (3.10) 

as as 

s=so 



which is negative when Pab is positive definite. This contradicts (|3.9|) and hence we 
conclude that our assumption was wrong and there are no closed timelike curves. 

For example, the stationary Godel type metric ( |2.1| ) represents a particular case of 
nj) with a = 1 and 
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^oo^ 



z/, = (0,v^,0), 



ab 



k 
1 



^ 0^ 



e"'^ 
1 



(3.11) 



It is clear that (3ab in ( |3.11|) is positive definite for k > 0. Negative values of k, as we 
saw explicitly in Sec. ^ allow for the closed timelike curves. In particular, this is true 
for the original Godel solution ( |2.9| ). 



Isotropy of MBR: The Microwave Background Radiation (MBR) is totally isotropic in 
( |3.lD for any moment of the cosmological time t. 

In the geometric optics approximation, light propagates along null geodesies, i.e. curves 
x^(A) with certain affine parameter A and tangent vector k'^ = dx'^/dX which satisfy 



rv^fc^ = 0, A;^ k. = 0. 



(3.12) 



As it is well known, see e.g. p5| , p6| , the red shift Z, which reflects the dependence of 
frequency on the relative motion of source and observer, is given by the simple formula 



l + Z 



{kf'u^,) 



(3.13) 



where the subscripts "S" and "F' refer, respectively, to the spacetime points where a 
source emitted a ray and an observer detected it. 



MBR is characterized by the black-body spectrum. One can show that the tem- 
perature of the black-body radiation depends on the red shift. 



1 



{kf'u 



(k'^u^) 



(3.14) 



In general, as was noticed, e.g., in ||142| , pT| , |32| , p3| , |25| , the observed temperature depends 
on the direction of observation through the ratio of (A;'^m^)'s. However, it is not the case 
for the class of shear- free models ( ^.1| ). Indeed, from ( |3.12[ ) we easily prove that the 
scalar product of the wave vector with any Killing vector or conformal Killing vector 
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is constant along the ray. In particular, for the conformal Killing (|3.8|) we immediately 
obtain the first integral of the nuU-geodesics equations ( ^.121) : 



fc/.Cnf = const =^ a{tp){u^k^)p = a{ts){u^'k^)s, (3.15) 

where we used the evident proportionality of the Killing vector to the average four- 
velocity of matter, 

Cnf = aWn^. (3.16) 



Substituting ( ^IsD into (133^ ), one finds 



Tp = -^ = T.i^M. (3,17) 

1 + Z a[tp) 



Thus, in the class of cosmological models (|3.1| ), MBR is completely isotropic, exactly 
like in the standard cosmology. As a consequence, the earlier strong restrictions of the 
global rotation derived from the isotropy of MBR |p.42| , ^ , ^ , ^ ,^ , are not applicable 



to the class of models under consideration. At the same time, the result obtained by 
no means contradicts the earlier observations, all of which were made for cosmologies 
with shear. Here we have trivial shear (|3.6|) , and thus one clearly concludes that it is 
shear, not vorticity, which is the true source of anisotropy of MBR. 



No parallaxes: Rotation (|3.5| ) does not produce parallax effects. 

In a very interesting paper [ |141|| , Hasse and Perlick introduced the generalized notion 
of parallax in cosmology. Let Si and 5*2 be any two sources, seen by an observer P. 
Provided all the three objects do not have peculiar motion with respect to the matter 
with the average four- velocity u^, one says of a parallax effect when the angle between 
the rays coming from 5*1 and 5*2 to P, changes in time. Six equivalent conditions were 



formulated in ||141|| which are sufficient for the absence of parallax. One can easily check 
that all of them are satisfied for the class of models ( |3.1| ). For example, according to the 
P(3) condition (in notation of ||141|| ), the velocity u^^ should be proportional to some 



conformal Killing vector. Indeed, this is the case, ( |3.16|) . Alternatively, according to 
P(4) there should exist a function / such that 

One verifies straightforwardly that / = Ina. According to P(5), is shear-free and 
the one-form 



P 



is closed, dp = 0. From ( |3.7| ) and ( |3.4| ), we find p = 0, hence P(5) is trivially fulfilled. 
Finally, according to P(6), there should exist a "red shift potential". One can prove 
that again the latter if given by the function / = Ina. 

We thus conclude, that the models ( ^.11 ) are parallax- free, and consequently, the value 
of vorticity cannot be determined from parallax effects, contrary to the claims of 
P73|J417|J4T8[ . 



Summarizing, cosmological models with rotation and expansion ( |3.1|) solve the first three 
problems of cosmic rotation, and the most strong limits on the cosmic rotation, obtained 



earlier from the study of MBR [142,B1,32] and of the parallaxes in rotating world [417,373 



are not true for this class of cosmologies. This class of metrics is rich enough, as it contains 
all kinds of worlds: open and closed ones with different topologies. 

4. GODEL TYPE EXPANDING COSMOLOGICAL MODEL 



To the end of this paper we will consider now the natural non-stationary generalization 



of the original Godel metric (|1 . 1|) which has drawn considerable attention in the literature. 
This generalized model is obtained from (2J^) by introducing the time-dependent scale factor 
a(t). 



ds^ = dt" - 2^a{t)e'^''dtdy - a\t){dx^ + ke^'^'^dy^ + dz^). 



(4.1) 
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Here, and the metric ( |4.1| ) evidently belongs to the family ( |3.1| ) with 

As we saw already, the condition A; > guarantees the absence of closed time-like 



curves. The metric ( [4 .11 ) is usually called the Gddel type model with rotation and expansion. 
Coordinate z gives the direction of the global rotation, the magnitude of which 



/I m a , , 

u; = W -cj.^cj^^ = — J- (4.2) 



decreases in expanding world [compare with ( |2.2| )]. 

The three Killing vector fields are, recall ( p.lO| ) and ( |2.16|) , 



^(1) = ^(2) = d^, ^(3) = —dx - ydy, . (4.3) 

These satisfy commutation relations 

[^(1), ^(2)] = ^(2), [^(1), ^(3)] = [^(2), ^(3)] = 0, (4.4) 

showing that the model ( |4.1| ) belongs to the Bianchi type III. 

As before, it is convenient to choose at any point of the spacetime ( [4.1| ) a local orthonor- 
mal (Lorentz) tetrad /i^. For simplicity, we will work with the immediate replacement of 



/iO = 1^ hO = -aV^e"^", h\ = hl = a, hi = ae"'''VkT^. (4.5) 

From the point of view of the Petrov classification, one can verify that the Godel type 
model is of the type D. Indeed, a straightforward calculation yields an obvious generalization 



of (|;6D: 



^Oioi ~ ^0202 ~ ^2323 ~ ^3131 ~ 2 ^i2i2 ~ 2 ^0303 ~ Qa'^(^t) + oj ' '^^'^^ 



5. DYNAMICAL REALIZATIONS 



Several dynamical realizations (i.e. construction of exact cosmological solutions for the 
gravitational field equations) of the shear-free models (|3.1|) , and in particular of the Godel 
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type metric ( [4.1|) , are known. In the Einstein's general relativity theory the Godel type 
models were described in |pi2| , |19(]| , pil| - |318| , p20|| , with different matter sources, whereas the 



Bianchi type 11 and 111 solutions (|3.1| ) were obtained in ||357| , |283|| for the imperfect fluid with 



heat flow and cosmological constant. 

In this section, we will describe a model [ p. 6 7] ] in which the gravitational field dynamics is 
determined by the minimal quadratic Poincare gauge model, which is the closest extension 
of the Einstein's general relativity theory. 

A. Field equations 

In the framework of Poincare gauge theory, gravitational field is described by the tetrad 
/i^ and the local Lorentz connection F";,^. In general case, gravitational Lagrangian is 
constructed as an invariant contraction from the curvature tensor 

R°'bfj.u = d^V^b^ — d^V^bfj. + ^'^c^i^'^bu — ^"'cu^'^bfi, (5.1) 

and the torsion tensor 

= d,K - + f\X - f (5.2) 
Note that here our notations for the connection, curvature and torsion are slightly different 



from [149]. In order to illustrate some common properties of cosmological models with 



rotation and expansion, here we consider the so called minimal quadratic Poincare gauge 



gravity [167] which is the closest extension of Einstein's general relativity (we may also call 
such a model the generalized Einstein-Cartan theory). The corresponding Lagrangian reads: 
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i^ + ^i^^ + aiTV^a'^' + aaT^^.T'^^^ + asT^T'^ . (5.3) 



Here R = h'^h'^^ R^biiu ^ the Riemann-Cartan curvature scalar, T°'^u = h'^T'^f^i,, and = 
T^^x - is the torsion trace, b, ai,a2, as are the coupling constants. 

Independent variation of the action ( p. 3] ) with respect to /i^ and P^f,^ yields the field 
equations: 
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+ a3 Q g^uTaT'' - V,T^ + g^uVaT''^ = (5.4) 

+ (2ai - aa) T[h" - a2T%, + ag^f^T,] = Kr%,, (5.5) 

where k = SttG, with r^^i^ and as the canonical tensors of spin and energy-momentum, 
respectively. 

Equation (|5.5| ) generalizes the well known algebraic relation between torsion and spin in 
the Einstein- Cart an theory. The tensors of torsion and spin may be decomposed into irre- 
ducible parts: trace and = T^f^x, pseudotrace = \e^vapT''°'^ and = \e^^apT^°'^ ^ 
and trace- free parts T"^,^ and t'^^u'- 

= + ^T[^5^] + ^e\,,T\ (5.6) 
T%. = + ^T[^6:^ + ^e%,,r\ (5.7) 

Easy to see that eq. ( |5.5| ) decomposes into three equations for irreducible components: 

(/i2 + 2bR) + 3bd^R = KT^, (5.8) 
(/ii - bR) = KT^, (5.9) 
(/i3-6i?)TV = 'tr%,, (5.10) 

where we denoted the following combinations of the coupling constants 

1 

/ii = 2ai - 2a2 - -, (5.11) 
/i2 = 1 - ^(2ai + a2 + 3a3), (5.12) 
/^3 = -Q + ai + |), (5.13) 
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which play important role in the determination of mass for quanta of torsion and couplings 
in all spin sectors of quadratic Poincare gauge gravity model. 

The eq. ( p.8| ) shows that in contrast to the other irreducible parts, the torsion trace is 
characterized by the differential and not merely algebraic coupling. The Riemann-Cartan 
curvature scalar plays a role of its "potential" . Contraction of ( |5.4| ) yields an equation for 
i?, 

-R+^{1 + 2/i3) T^^.T"'^' + ^ + ^/^^"^ + 2(1- /X2) (V^T'^ - ^T^T^) = kE, 

(5.14) 



where S = g^^pTj^^'^ is the trace of the energy- momentum tensor. Inserting ( |5.8|) -( pn0 ) into 



( p.l4| ), we find a non-linear differential equation for the curvature scalar. As usual, hereafter 
we denote by tilde the purely Riemannian (torsion-free) geometrical objects and operators. 

In view of specific nature of the torsion trace, it is convenient to separate its contribution 
to the eq. (^^). Let us introduce 

T"" = T + -^'^fiuxT^, (5.15) 

and hereafter hat will denote the corresponding geometrical quantities constructed with 
the help of the Riemann-Cartan connection with the trace- free torsion ( [5.15| ). The symmetric 
part of ( |5.4| ) then is rewritten as 

(1 + 2bR) [r^^,,) - ]^g,,R) + hn^g,, + (^f^) T^%.)x 

+ ( ^^^3 ^^' ) (V(^T.) - g,.VaT^ - \g,.T^T^ ) = (5.16) 



whereas the antisymmetric part of ( ^.4|) is an identity in view of the geometric relation 

2i?[H + (V„ + (t%, + 25f^T,]) = 0, (5.17) 
and the angular momentum conservation law 
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(V„ + T„)rV = S[H- (5.18) 



Inserting the torsion and T from ( |5.9| ) and ( p.lQ ) into ( p.l5| ), and further into ( p. 161 ) 



one can write, similarly to the Einstein-Cartan theory, the effective Einstein equations for 
the metric. Spin of matter will contribute to the so called modified energy-momentum 
tensor. 

B. Matter sources 

We will describe cosmological matter by means of phenomenological model of a spin- 
ning fluid of Weyssenhoff and Raabe. This is a classical model of a continuous medium 
the elements of which are characterized, along with energy and momentum, by an intrinsic 
angular momentum (spin). Usually, hydrodynamical approach is considered as a good ap- 
proximation to the description of realistic cosmological matter on the early, as well as on 
the later stages of universe's evolution. It seems natural to take into account spin (proper 
angular momentum) of elements of cosmological fluid, i.e., particles on early stage and galax- 
ies and clusters of galaxies on later stages. In general case, an element of the Weyssenhoff 
fluid is characterized by the tensor of spin density r^^,, charge e and (proportional to spin) 
magnetic moment xr^jy (x - is a constant). The consistent variational theory of spin fluid 



in the Riemann-Cartan spacetime, developed in ||289|| , yields the following currents in the 



right-hand sides of the gravitational field equations (|5.4| ) and 

-F^^F,'' + ^g^^F^pF^^ - 2x (r^.r," + r^pF^^u^u^) , (5.19) 
= m"v, (5.20) 

where is fluid's four- velocity, p is pressure, e is the internal energy density, F^^, = 
dfj,Ai, — duAfj_ is the electromagnetic field created by and interacting with the electrody- 
namical characteristics of the medium. 

Dynamics of the fluid is given by the rotational and translational equations of motion: 
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[{S^^ - uf'u,) r,, - (5^ - u^u,) , (5.21) 
(V. + T,) S% - T^^S'^, - T^pR'^^^, = 0. (5.22) 



The self-consistent variational principle ||289|| yields also the Maxwell's equations for the 



electromagnetic field and the conservation law of number of particles: 

V^F^" = epu" + 2xV^r^^ (5.23) 
{pu^) = 0, (5.24) 

where p is the particle density. 

The system of differential and algebraic equations { WM i ( |5.5D , ( |5. 191) -( ^^241) is complete, 
and one can determine all the gravitational and material physical variables from this system. 
It is worthwhile to remind that spin satisfies the Frenkel condition, 

u^T^, = 0. (5.25) 

The specific features of dynamics of the Weyssenhoff spin fluid in the Riemann-Cartan 
spacetime were studied in detail in ||289 |. 



C. Exact solution for the Godel type universe 

Let us describe cosmological model with rotation and expansion which provides an exact 
solution of the gravitational field equations. We will restrict ourselves to the shear-free class 
of spatially homogeneous metrics ( p.l| ), and for definiteness, a particular case - the Godel 
type model ( [4.1|) - will be considered in detail. 

The field equations ( p.4| ) and ( |5.5| ) with the sources ( p.l9| ), ( p.20| ) are most conveniently 



formulated with respect to (nonholonomic) local Lorentz frame. Let us choose the coframe 
( [4.51 ). Its inverse reads: 



hi = 1, hl = Jr^, hl = hl = -, hi = \^=, (5.26) 
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Direct calculation of the local Lorentz connection F^h^ = h'^h'^Tp^ + h'^d^h'^ (with Tp^ as 



the Christoffel symbols) yields for the metric (jO 



20 ~ 21 — 23 ~ ^ y T , ' ^ 22 — ^ ■> ) 

a V K + a a 



ni-r^o--r°i2-£y^, (5.28) 

where f\, = T\^h^ . 

Since the direction of rotation (along the z axis) is apparently a distinguished one, it 
is natural to assume that the spin of the fluid and electromagnetic field are also oriented 
correspondingly along z. For the medium with vanishing magnetic moment of its elements 
(X = 0) one finds from (p:2ll)-(p:25|): 

ri2 = roa-3, (5.29) 
F^2 = fioa~^ (5.30) 
p = Poa"^ (5.31) 

where the integration constants are related by the consistency condition 



mBo J ^ =epo. (5.32) 
\j k + a 

We will now describe a solution of the gravitational field equations for which the Riemann- 
Cartan curvature scalar is assumed to be constant: 

R = const. (5.33) 



Then in view of ( |5.25| ), eqs. ( p.8D -( pTTUD yield the irreducible parts of torsion: 



= 0, T^ = ye^,,^MV°^ T"^, = ^ (2m-v + - «-^%) . (5-34) 
where we denoted the constants 

Ai = A3 = (5.35) 

Pi — oR Pa — oR 
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Nontrivial solution exists only when /ii ^ hR ^ /is, otherwise ( |5.9| ) and ( |5.10|) lead to the 
vanishing spin. As for the eq. (|5.8| ), one has two alternatives: if ^2 + 26-R 7^ 0, then the 
torsion trace is zero, whereas if fi2 + 2bR = 0, then is arbitrary. However in both cases the 
torsion trace completely decouples from the eq. ( [5.16| ), hence without restricting generality 
we can put = 0. 

Using (|5.34|) , we can rewrite (|5.15| ) as 
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= - [(2A3 + Ai) M°r^, + (A3 - Ai) M^r% + (Ai - A3) M,r%] . 



(5.36) 



Substituting this expression into ( ^.16| ), one finds 



= (2A3 + Ai) [Ai (1 + 2/ii) - A3 (1 + 2/^3)] 

io 



+ 



'Ai -4A3 
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[Ai (1 + 2/ii) + 2A3 (1 + 2/i3)] u^u.r^pr'^^ 



+ 
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4A^ (1 + 2/i3) - A^ (1 + 2/ii) g^, r^pr'^^ + A3 (1 + 2/13) V„ u^^t,) 



(5.37) 



After these computations, one finally finds that when the constant in (|5.33|) is equal R = 
—1/26, the equations ( ^.16| ) are reduced to the algebraic relations between the pressure p, 
energy density e, spin Tq,^ and Maxwell tensor Fq,/3. From these equations, in addition to 
( ^.32| ) we find a new relation between the integration constants. 



Bl = -niTQ 



k + a 



(5.38) 



Without restricting generality, one can choose the parameter m to be positive, then the 
minus sign in (|5.38| ) means that the spin and vorticity have the same direction. Thus tq < 



[from dX^), and (§;^ we have 
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la v fc+^]- turn, the direction of the 



magnetic field (sign of Bq) depends on the sign of the charge e. 

With the help of ( ^.38|) the equations ( p.l6| ) are reduced to the following two equations 
which describe the state of matter: 
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8bK ' 2a^ ' \ 6 y a6' ^^'^^"^ 



Finally, the evolution of the scale factor a(t) is given by the equation (|5.33|) . Substituting 
(PD, ( p:26D , and ( p:28D in it, one finds 



a a 



2 



As one can see, depending on the values of the coupling constants A3 and Ai, torsion can 
either accelerate or prevent the cosmological collapse. We will consider the latter possibility, 
assuming, for concreteness, that 4A3 ^ Af. 
Integrating ( ^.41| ), one gets 



^2 ^ m\3a + Ak) a\k + a) _ {k + a) (4Ai - A^) Tq^ a 

12k 2Akb 36ka^ a^' ^ ' ' 

The integration constant a is determined by the normalization of the scale factor a = 1 at 
the moment of a bounce (t = 0), when a = 0, 



a 



2 ^2^ 2 (3cr + 4fc) k + a 



D. Estimates of parameters of the cosmological model 



In order to integrate the equation ( |5.42| ), one needs to know the constant parameters 



which enter that equation, because their numeric essentially affect the form of solution. All 
constants are uniquely determined by the data which describe modern state of the universe: 
the equation of state ~ (dust), energy density ct, Hubble constant Ht{= a/a) and the 
deceleration parameter qt (= —aa/a?) are estimated by 

0.15 < ^< 12, 
4- 10"^ Vr"^ < Ht< lO^^Vr^S 
0.01 < qt<l. 
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(Hereafter the index t denotes modern values of the physical quantities). 

Previously, the angular velocity was estimated as relatively small oJt/Ht 10^'^. 
New analyses of observational data suggests higher estimates for the global vorticity, more 
close to Ht- For the qualitative analysis of the model let us choose 



qt = 0.01, = 0.18, = 0.1. 



(5.44) 



Besides the estimates (|5.44|) , one needs also the values of the spin-torsion coupling constants. 
There is a large ambiguity which affects the scale of the space at the moment of a bounce. If 
one assumes, cf. Trautman |[416|] , that the cosmological collapse stops on the scale of order 
of 1 cm, then we get from ( ^.44 ) an estimate 



K 



2.3-10 
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(5.45) 



Ai - 4A3 

Comparison of the equations ( p.41| )-( p.43| ) with ( p.44| )-( p.45| ) yields the modern value of the 
spin density 



3 ■ lO^gcm-^-^ < n < 7.6 ■ 10^gcm-^s-\ 



(5.46) 



the geometric parameters of the model 



-^71, 
a 



(5.47) 



the coupling constant 



(5.48) 



and the magnetic field 



10-^G <Bf< 2.6 • 10-^G. 



(5.49) 



The final integration of the equation ( ^.42|) determines the law of evolution of the scale 
factor: 

2/i 



A^/a^Htt= (const) + InZ(a^) + 
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7i 
7i + 72 



F{(^] m) — CiU.{lp] n, m) 



(5.50) 



where the integration constant (const) is such that a = 1 for t = 0, and we denote the 
function of the scale factor 

(7i+72)(7i-72 + 2a2) 2 
In eq. ( ^.50|) , the constant parameters are determined by: 



m^{?>a + Ak) k + a\ AXl-\\\ k + a 



«i = TTTTTTp ' "2 = — ; — — 779' "3 



12kH! ' ' \ k J \ 36 J nr l2hkHr 
ui > U2 are the roots of the quadratic equation 



+ [1 + — ]u + — =0, 
as J as 



from which one constructs 



A± = -U1U2 (^y 1 - i/ui ± - 
_ A+(l-A-) _ A_ (A+ - 1) _ 1-A„ ^ _ A+ - 1 

-Dl — ^ ^ , Oi — r r , J32 — 1 ^ — 7 ^2 



X-L^ — A_ A_L — A_ A^ — A_ A_L — A_ 



/ U1U2 I U1U2 2 Bi A+ 

= (A+ - A_) /A+, n = Bin? = 1 - A^. 

Finally, F {if; m) and Il{ip; n, m) are elliptic integrals of the first and the third kinds, respec- 
tively, with their argument given by 

/ 2 \ 2 

2 1 -2 / « - 72 \ 



(l — k sin^ ip = 1 — 1^ 



a? - 7i 



The complete picture of the cosmological evolution is obtained when one substitutes ( |5.44|) - 
( p.49| ) into the above equations. In particular, one can verify that the age of the universe 
[estimated as the time since t = till the modern stage characterized by ( |5.44| )] is equal 
T ^ H^^^ in a good agreement with observational data. 

The analysis of the evolution of the scale factor a(t), see ( [5.50| ) and ( [5.42| ), shows that one 
can naturally distinguish several qualitatively different stages in the history of the universe. 
First stage is the shortest one and it describes a bounce in the vicinity of t = 0. There is no 
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initial singularity due to the dominating spin contribution in ( p.42| ). At this stage, the fluid 
source is characterized by the approximate equation of state of stiff matter 



P 



' Ai-4A3 \ 
6 / a6 



The applicability of the classical (non-quantum) gravitational theory is guaranteed by the 
condition on the curvature invariant Rai3nvR"^^^ ^ ^p/, which is satisfied at the moment of 
the bounce. The duration of that stage is quite small (-C 1 s), since the spin term quickly 
decreases in (^.42|) with the growth of the scale factor. 

Next comes the second stage when the scale factor increases like v^, while the equation 
of state is of the radiation type, p ~ £:/3. This "hot universe" expansion lasts until the 
size of the metagalaxy approaches ~ 10^"^ cm. After this the "modern" stage starts with 
the effectively dust equation of state pt ~ 0, ~ Scale factor still increases, 

but the deceleration of expansion takes place. The final stage depends on the value of 
the cosmological term, and either the future evolution enters the eternal de Sitter type 
expansion, or expansion ends and a contraction phase starts. 

Our purpose in this section was to demonstrate the exact solution which avoids the 
principal difficulties of old cosmological models with rotation and expansion. Certainly, one 
cannot claim that the results obtained describe the real universe, and consequently we have 
limited ourselves to a test type computations for the physical and geometrical parameters, 
without trying to find the best estimates. The main difficulty of this particular dynamical 
realization, in our opinion, is presented by the magnitude of magnetic field (|5.49|) which at 
the "modern" stage should be close to the upper limits established for the global magnetic 



field from astrophysical observations, see e.g. |24 



It seems worthwhile to note that ( [5.46| ) gives a rather big value for the spin density 
10^^^/cm^). This is in a good agreement with the macroscopic interpretation of the 
cosmological spin fluid, the elements of which are spinning galaxies. Indeed, Tt/ pt ~ 10^°^/l 
is close to the observed magnitude of the angular momentum of a typical galaxy. Moreover, 
for the total spin of the metagalaxy occupying the volume V ~ lO^^cm'^, one finds rtot = 
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TtV ~ 10^^°^. This once again draws attention to the Large Number hypothesis of Dirac- 
Eddington and related questions, see [p'^- pO|j63|J44CI|| and ||271| - |274|| . 



6. NULL GEODESICS IN THE GODEL TYPE MODEL 

Practically all the information about the structure of the universe and about the prop- 
erties of astrophysical objects is obtained by an observer in the form of different kinds of 
electromagnetic radiation. Thus, in order to be able to make theoretical predictions and 
compare them with observations, it is necessary to know the structure of null geodesies in 
the cosmological model with rotation. All the models from the class have three Killing 



vectors and one conformal Killing vector field. Hence the null geodesies equations (|3.12|) 
have four first integrals, 

% = Cnf^M' (la= - ^(1)^^' a = 1, 2, 3. (6.1) 

Solving ( |6.1D with respect to k^, one obtains a system of ordinary first order nonlinear 
equations which can be straightforwardly integrated. Complete solution of the null geodesies 
equations in the Godel type model is given in ||186| , P87|| , and here we present only short 
description of null geodesies in ( ^4.1|) . 

To begin with, let us define convenient parametrization of null geodesies. Without loosing 
generality (using the spatial homogeneity) we assume that an observer is located at the space- 
time point P = {t = to,x = 0,y = 0, z = 0). Now, arbitrary geodesies which passes through 
P is naturally determined by its initial direction in the local Lorentz frame of observer at 
this point. In the tetrad ( |4.5| ) we may put 

k% = {hlk^')p = (1, sin ^ cos 0, sin ^ sin 0, cos^), (6.2) 

where 6, are standard spherical angles parametrizing the celestial sphere of an observer. 
Then from ( p.8| ), ( [4.3| ), ( |6.2| ) and ( |6.1| ) one finds the values of the integration constants 

go = ao, 
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Qi = CLoiV^ + V/c + cr sin ^ sin ( 



(6.3) 



q2 = ao cos 9 
ao 



93 



sin 9 cos (b. 



m 



Generic null geodesies, initial directions of which satisfy gi/ao = (sin 6' sine 
are then described by 



+ 



fc+o-. 



y/a + y/k + a sin 9 sin $ 
y/a + \/k^ G sin 6' sin ' 
sin ^( cos $ — COS0) 

771(^/0"+ + crsin6'sin< 

t 



k + a 



cos 9 



lio 



dt' 



+ 



a 



k + a 



m 



where the function $(t) satisfies the differential equation 



'dt 



m 
a 



■jrr- + sin 9 sin $ 



1 + Jfc^sin^sin<l> 



(6.4) 

7^0, 

(6.5) 
(6.6) 

(6.7) 



(6.8) 



with initial condition $(to) = 0- 

For a detailed discussion of rays which lie on the initial cone (sin 9 sin + , / = and 
different subcases of (16. 51) -(16. 81) see [E87 . 



7. OBSERVATIONS IN ROTATING COSMOLOGIES 

The qualitative picture of specific rotational effects which could be observed in the Godel 
type model ( [4.1|) is in fact independent of the dynamical behaviour of the scale factor a(t). 
To some extent the same is true also for the quantitative estimates, especially if one uses 



the Kristian- Sachs formalism ||2 14| , |240|| in which all the physical and geometrical observable 
quantities are expressed in terms of power series in the affine parameter s or the red shift 
Z . In this case the description of observable effects on not too large (although cosmological) 
scales involves only the modern values (i.e. calculated at the moment of observation t = to) 
of the scale factor oq = a(to), Hubble parameter Hq = {d/a)p, rotation value uq = uito), 
deceleration parameter go = — (aa^/a^)p, etc. 
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In this section we will outline possible observational manifestations of the cosmic rotation 
in the Godel type universe. Estimates for the value of vorticity and for the direction of 
rotation axis can be found from the recent astrophysical data, see Sects. [7 C| and |7D. 



A. Kristian-Sachs expansions in observational cosmology 

Kristian and Sachs in the fundamental paper |pi4|| have developed universal technique 
for the analysis of the observational effects in an arbitrary cosmological model. Essentially, 
their method is based on the theory of null geodesic (ray) bundles in Riemannian spacetime. 
The main advantage of the Kristian-Sachs approach is in its generality which makes it 
possible to consider cosmological tests even without the knowledge of explicit solution of 
the gravitational field equations. This is particularly convenient when one is interested 
in revealing additional (as compared to the standard Hubble expansion) geometrical and 
physical properties of the universe. In this section we will apply the method of Kristian and 
Sachs to the rotating cosmological models. 

Let k^^ be a wave vector field which satisfies the null geodesies equation ( |3.12| ) in a 
Riemannian manifold. (The connection is torsion-free, and we will omit the tilde to simplify 
the notation). The starting point is the equation 

= -R'^^^pA'^k'^k^, (7.1) 

which determines the deviation vector A'^ := 6x'^ of the rays in a bundle. Hereafter dots 
denote covariant derivatives along a ray: (') = k'^'Vfj,, e.g., A'^ = DA^/ds = k^'V^A'^, etc. 

The crucial idea of Kristian and Sachs is to find a solution of ([7.1|) in the form of a 
series in the affine parameter s. Correspondingly, all the physical and geometrical quantities 
are ultimately represented by series in s. Let us give some basic results of that approach, 
without going into calculational details. 

We will denote the spacetime points of an observer and a source by P and S, respectively. 
Let us start with the red shift defined by ( ^.13| ). A simple Taylor expansion of the right-hand 
side yields: 
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7 ^ (3) 

2 6 



Here we denoted: 



-s := s {k^u,)p, = (^-^^ , := {k^^ . . . K^-'^K^"V,, . . . V,„_,«,„^^ 

Note that the value s describes a position of the observer on a null geodesic connecting S 
and P, so that s = corresponds to S. The minus signs in ( [7 .21 ) arise when we expand 
{u^'k^)s at P. 

Let us consider a bundle of rays emitted from a point source 5*. We are interested in 
A^(s) along a geodesic from 5* to P. One has A'^(O), and hence from ( [7.1D , A^(0) = 0. 
The vector A'^(O) is chosen to be orthogonal to k'^ and u'^ at S, i.e. we are studying the 
cross-section of the bundle which is orthogonal to the central ray in the rest frame of the 
source. Note that such an initial condition provides orthogonality k'^A^ = at any value of 



s, in view of (fO 



The area distance is one of the central notions in the observational cosmology. If an 
observer P is viewing the distant source S with an intrinsic perpendicular area dAs which 
subtends the solid angle dflp at P, the area distance r between P and S is defined by 

dAs = r^dnp, (7.3) 

Another (reciprocal) area distance between S and P is defined by dAp = r'^dQs^ where 
the solid angle dfls subtends the area dAp at P. The reciprocity theorem Il95|j96| , [240| tells 
that these distances are related by 

r| = (1 + Zf. (7.4) 

One can obtain vs with the help of A'^(s) which describes the behaviour of the bundle of 
rays. 

The deviation equation ( |7.1| ) determines iteratively the third and higher covariant deriva- 
tives of A'^ at S, and thus one construct the solution in the form of the Taylor series 

A^(s) = sA^(0) + ^A^(0) + ... . 
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In order to find r^, one needs actually only the length of A'^. Straightforward calculation 
yields: 



3 _ 



A^(s)A^(s) = s2A'^(0)A^(0) 

+ A'^(0)A'^(0) 
The spacetime curvature contributes via the expressions 

h '.= R^u k , Cq^ := Cdfif^u k^ k 
higher order corrections are contained in the terms 0{s'^). 



J 5 



(7.5) 



(7.6) 



After some algebra, one finds from (|7.5| ): 



6 12 ^ ' 



(7.7) 



J 5 



Combining ( ^.1^ ) , (|7^) , and ( |7.7| ) , one can easily derive the series which represents the area 
distance r in terms of s. 

However, since the value of the affine parameter s for the source S and the observer P 
cannot be measured directly, the series in powers of s are not particularly useful. Instead, 
Kristian and Sachs proposed to invert the relation ( [7.2|) , and then substitute s = s{Z) in all 
other quantities. Direct calculation yields the following final series representation: 



l + Z= l- ru + Y u - — lu+-Bj+..., 

where, replacing ([7. 61), B := {R^^ K'^)p. Inverting ([7.8D, one finds 
Z 



(7.8) 



(1) 
— u 



z 


(2) 
U 


Z2 


r (2)2 


1 u 


2" 


u ^ 






~ 6 (1)^ 



u ■ 



(7.9) 



B. Classical cosmological tests 



Classical cosmological tests, such as apparent magnitude - red shift (m — Z\ number 
counts - red shift (A^ — Z\ angular size - red shift relations, and some other, reveal specific 
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dependence of astrophysical observables on the angular coordinates {0, 0) in a rotating world. 
Thus a careful analysis of the angular variations of empirical data over the whole celestial 
sphere is necessary. 

The knowledge of null geodesies enables one to obtain the explicit form of the area 
distance r between an observer at a point P and any source S. Each ray (a null geodesic) 
x^^s-jP"^) passing through P is labelled by the two parameters := 0), A = 1,2 which 
have a simple meaning of the spherical angles fixing the direction of ray at P. Thus in 
general, r is a function of the direction of observation, that is r = r{6, 0). Besides this, r of 
course depends on the value of the affine parameter s, or equivalently on the moment of t 
at which source radiates a ray detected by an observer at to. The area of the cross-section 
of the bundle of rays is 

dx^ dx^ 

QpA Qpl 



dAs = \/det Gab d6d(j), Gab = 9i,u tt^- (7.10) 



Since the solid angle in ( |7.3| ) is defined, as usually, by dQp = sin 9 dOdcp, one easily obtains 
the general expression for the area distance between an observer P and a source S which 
lies on the geodesic labeled by the observational angles {9, 0) at the value s of the geodetic 
parameter. 



Jdet GAB{s\9,(t)) 
sm 9 

The knowledge of the null geodesies (|0|)-(|0|) makes it possible to calculate area distances 
in the Godel type universe explicitly. For example, using ( |6.5| )-( |6^ ), one can find for the 
area distance along the axis of rotation the following exact result 

sin^ (jdt'uj{t')\ 
At;9 = 0) = ^. (7.12) 

UJ^t) 

The rotational effects are maximal in the directions close to the z axis, and quite remarkably 
the observations in the direction of rotation can be described by simple and clear formulas. 

As for an arbitrary direction, exact formulas become very complicated and it is much 
more convenient to replace them by the Kristian-Sachs expansions. One can use the ex- 
pansions (|7.8| ) and ( [7.9|) in the calculations of observable effects in rotating cosmologies: 
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All the angular dependent rotational contributions are contained in the coefficients of these 
expansions. 

For the Godel type cosmology ([4.1|) the classical {m — Z) and (A^ — Z) relations read as 



follows p7| , p88|p7| 



Apparent magnitude m versus red shift Z: The apparent magnitude of a source 5* 
with the energy flux Lp, measured at P, is defined by m = — |log]^o ^p- Given the 
total output of the source L, Ls = L/Att, the luminosity distance is := Ls/Lp = 
r^(l + Z)'^, see [p5ip6| . Using ( [7. 91) , one finds for the Godel type universe: 



m = M -5 logio ^0 + 5 logio Z + - (log^o e){l-qo)Z 



-Slogio 1 + 



k + a 



sin 6 sin i 



- ? (logio ^ ^ sin 9 cos 



^ + sin^sin, 



7-2— sin 9 sin i 

k+<T 



0{Z^). 



(7.13) 



Here M = — llog^^Q-L^ is the absolute magnitude of a source with an intrinsic lumi- 
nosity Ls- 

Number of sources A^ versus red shift Z: Let us assume that the sources of electro- 
magnetic radiation are distributed homogeneously with the average density n(t). Ne- 
glecting the evolution of sources, one can count the number of images seen by an 
observer as follows. Consider the bundle with a solid angle dQ detected by an observer 
at P. When moving along the central ray, the distance between the two infinitesimally 
close cross-sections is equal di = {k^u^)ds, and the volume between these sections is 
dAd^, with dA = r'^dQ. Hence, the number of sources contained between these two 
sections is given by dN = —ns{k'^u^)sr^dflds. Integrating, and using ( [r.8|) -( [r!9D : 



dN 



n^Z^ 



3HI [1 + Jj^ sin sin, 



'0 V " ' \l k+a 

sin 9 cos d) 



l--{l + q,)Z 



-3 



sin 9 sin ( 







(l + sin sine 



Z + 0{Z'^) 



(7.14) 
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The (A^ — Z) relation ( |7.14| ) describes the number of sources per sohd angle dfl observed 



up to the value Z of red shift. One can estimate the global difference of the number of 
sources visible in two hemispheres of the sky, N^,N', by integrating ( |7.14| ). The result is 



- A^- 1 








k + a 


k + a) 



OiZ^). (7.15) 

It seems worthwhile to draw attention to the absence of a correction proportional to Z in 
( [r.l5| ). It is difficult to make a comparison of our results with [[431|J265| , |103|| as they study 



stationary rotating models in which there is no red shift. 

For some time classical cosmological tests were carefully carried out for standard models, 
but later it was recognized that evolution of physical properties of sources often dominates 
over geometrical effects. However, specific angular irregularities predicted in rotating cos- 
mologies, ( |7.13|) , ( |7.14| )- ([7T5|) , may revive the importance of the classical tests. 



C. Periodic structure of the universe 



Recent analysis of the large-scale distribution of galaxies has revealed an apparently 



periodic structure of the number of sources as a function of the red shift. Cosmic rotation 
may give a natural explanation of this fact ||193|| . The crucial point is the helicoidal behaviour 



of the null geodesies ( |6.5D -( p77D in the Godel type model in the directions close to the rotation 
axis. This yields a periodicity of the area distance as a function of red shift, and hence the 
visible distribution of sources turns out to be also approximately periodical in Z. 

This effect is most transparent for the direction of rays straight along the axis of rotation. 
The area distance is then given by ( |7.12|) . In order to be able to make some quantitative 
estimates, let us assume the polynomial law for the scale factor. 



a{t) = ao(- ^) , (7.16) 

\tO — too/ 

which is naturally arising in a number of cosmological scenarios (0 < 7 < 1). Then one can 
derive (analogously to ( [7.14D ) the distribution of number of sources per red shift per solid 
angle: 
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dN riQ , 2 ( l^Q 



sm 



(1 + - ll . (7.17) 



dndZ ulH^il + Zyh \{l-j)Ho 
Thus, the apparent distribution of visible sources is an oscillating function of red shift, with 
slowly decreasing amplitude. A similar generalized formula can be obtained for arbitrary 
directions, so that ( [7.17| ) is modified by additional angular dependence of the magnitude of 
successive extrema of distribution function. 

The observational data [HI give for the distance between maxima the value 128/i"^ Mps 
(where Hq = 100/ikm sec~^Mps~^). From this one can estimate the rotation velocity which 
is necessary to produce such a periodicity effect, 

ujq^IAHq. (7.18) 

This result does not depend on 7. 



D. Polarization rotation effect 

The cosmic rotation affects a polarization of radiation which propagates in the curved 
spacetime ( ^4.1|) , and this produces a new observable effect which has been reported in the 
literature by Birch |46.4^ and more recently by Nodland and Ralston | 279 , p80|| . In the 
geometrical optics approximation, polarization is described by a space-like vector which 
is orthogonal to the wave vector, /^fc^ = 0, and is parallelly transported along the light ray, 



(7.19) 



The influence of the gravitational field of a compact rotating massive body on the po- 
larization of light was investigated by many authors. 
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18,20,53, 


101, 


115 


,171 


258 



cosmologies shows that the global cosmic rotation also forces a polarization vector to change 
its orientation during the propagation along a null geodesies ||287| , |257| , |216|| . It is clear that 
this conclusion has physical meaning only when one can define a frame at any point of the 
ray with respect to which polarization rotates. Let us describe how this can be achieved. 
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As it is well known, gravitational field affects all the properties of an image of a source, 
such as shape, size and orientation | p77|j286|j95| , |96|j389[] . Like the rotation of the polarization 
vector, the deformation and rotation of an image depend both on the local coordinates and 
on the choice of an observer's frame of reference. However, if one considers a combination 
of the two problems, then a truly observable effect arises which is coordinate and frame 
independent. Putting it in another way, one should calculate the influence of the cosmic 
rotation on the relative angle rj between the polarization vector and the direction of a major 
axis of an image PI . This probl em was discussed in the papers [pl4|J141|| , but in our 
opinion, their results are incomplete, in the sense that only null congruences converging at 
the point of observation were considered. 

The relative angle r] can be most conveniently defined within the framework of the 
Newman-Penrose spin coefficient formalism. Namely, it is sufficient to construct a null 
frame {l,n,m,rn} is such a way that / coincides with the wave vector k, and the rest of 
the vectors are covariantly constant along /. Then we can consider m as a polarization 
vector, and thus a deformation of an image of a source, calculated with respect to the frame 
m, m}, gives the observable relative angle 1]. Let us describe explicitly the null frame: 
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(7.20) 
(7.21) 



(7.22) 



(7.23) 



Here 
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^(t, z) = z 



m 



a 



(7.24) 



2 V A; + (T 

One should note that (|7:20D -(|7:23D is a smooth field of frames which covers the whole space- 
time manifold. Direct computation proves that ( 7.20| ) is the null geodesic congruence with 
an affine parametrization. We say that this congruence is oriented along a direction given 
by the spherical angles {6, 0) in the local Lorentz frame of an observer at P, because a null 
geodesies (|6.5|)-(|0|) belongs to this congruence. 

It is straightforward to find the spin coefficients (we are using the definitions of [71], and 
mark the spin coefficients with tildes in order to distinguish them from the other quantities 
in our paper): 



£ = 0, 
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cos 6^ - 2 cos ^fcos 9-1] 



(7.25) 



(7.26) 



(7.27) 



We do not write the other spin coefficients, because their values are irrelevant. Only one 
important step remains to be done: the spin coefficient vf 7^ in the frame ( [7.20| )-( [7.23| ), 
and we need to make an additional Lorentz transformation 



/, n 



n + A*m + Am + A* A L m — > m + AL m 



m + A*L 



(7.28) 



where the function A satisfies equation l^dfj^A 



TC 



0. This ensures that vr = in a 



new frame, while remarkably the transformation ( |7.28D does not change any of the spin 
coefficients ( |7.25|) - (|7.27|) . Thus one obtains finally the field of null frames {l,n,m,m} with 
the required properties: / is the null geodesic congruence with affine parametrization, while 
n, m, m are covariantly constant along /. The latter is equivalent to k = e = vf = 0. 
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As it is well known, deformation and rotation of an image along a null geodesies are 



described by the optical scalars [p77 



= — Rep, u = Imp, 



a. 



(7.29) 



We now assume, for definiteness, that the polarization vector /'^ coincides with the vector 
of the above null frame. Let an image of a source, as seen at the point corresponding 
to a value s = si of the affine parameter, be an ellipse with the major axis p and the minor 
axis q. Then one can straightforwardly obtain [|194|| for the angle of rotation of the major 
axis of the image at S2 = si + 6s, 



_ p + q _ 

6rj = — ljj6s Imo" 6s. 

pi — q^ 



(7.30) 



Integration along a ray gives the finite angle of rotation. 

It is worthwhile to notice that exactly along the cosmic rotation axis the observer at P 
finds for the optical scalars 



^9/ 



0, 



Up — UJq, 



a = 0. 



(7.31) 



Thus the effect of rotation of the polarization vector in this direction is the most explicit. 
As for an arbitrary direction of observation, we finally find from (|7.30|) for k/a with 



the help of the Kristian-Sachs expansions (|7.71 )-([7 



T] = ujQr cos 9 + 0{Z'^). 



(7.32) 



This result is in a good agreement with the observational data reported |^6[ on the dipole 
anisotropy of distribution of the difference between the position angles of elongation (the 
major axis) and polarization in a sample of 3CR radio sources. The estimate for the direction 



and the magnitude of the cosmic rotation, obtained in p88| , |187| , |194|] from Birch's data, reads 



r = 295° ± 25°, 6° = 24° ± 20° 
cuo = (1.8±0.8)i/o- 



(7.33) 
(7.34) 
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Not entering into the discussion ||428|,|9^,p29|,|68| of the statistical significance of the results 



of Nodland and Ralston, we can make a comparison with the new data. In [ p79|| , the dipole 
effect of the rotation of the plane of polarization was reported in the form (3 = \A^^ r cos 7, 
where (3 is the residual rotation angle between the polarization vector and the direction of the 
major axis of a source (remaining after the Faraday rotation is extracted), r is the distance 
to the source, 7 is the angle between the direction of the wave propagation and the constant 
vector s. The analysis of the data for 160 radio sources yielded [p79|] the best fit for the 
constant = (1.1 ± 0.08) | m yr~^ (with the Hubble constant H), and the s'-direction 
RA(21h±2h), dec (0° ± 20°). In an attempt to explain their observations, Nodland and 
Ralston concluded that it is impossible to understand such an effect within the conventional 
physics. Instead, they considered a modified electrodynamics with the Chern-Simons type 
term violating the Lorentz invariance [0, and related to the coupling constant of that 



term. However, interpreting the effect as arising from the cosmic rotation, we immediately 
find 

/° = 50° ± 20°, 6° = -30° ± 25°, (7.35) 
Wo = (6.5 ± 0.5) ifo- (7.36) 

This is larger than the estimate (|7.34|) obtained from Birch's data Also the direction 



of the cosmic vorticity is different from ( [7.33| ). However, it is interesting to note that within 



the error limits, the two directions are orthogonal to each other. 

It is clear that further careful observations and statistical analyzes will be extremely 
important in establishing the true value of the cosmic rotation. In particular, it has been 



recently claimed in ||6^ , that the magnitude of the polarization rotation effect should be one 



or two orders lower than reported in ||279|| . The value of the cosmic vorticity then should be 
reduced correspondingly. 
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8. CONCLUSIONS 



In our discussion of the properties of rotating cosmologies we have paid special attention 
to the Godel type model ( |4.1| ). However the main conclusions are true also for the whole 
class of metrics One can notice, that in some cases the numerical estimates obtained 

for the value of the vorticity do not agree with each other, e.g., see ( |7.18| ), ( |7.34| ), and ( [7.36[ ). 
Such inconsistencies can reflect the fact that too few empirical data were analyzed until 
now, so that a further detailed discussion is required in order to obtain the final convincing 
estimates. Already now certain modifications of the above simple models are necessary: 
It is clear, in the light of the modern COBE results ||402| , |4l|j436[| , that the purely rotating 



models (|3.1|) should be replaced by the cosmologies with a nontrivial shear. A preliminary 
analysis of such generalizations shows that rotating models can be made compatible with 
the COBE data without destroying the rest of the rotational effects (in particular, without 
essential modification of the polarization rotation formulas). It may turn out, though, 
that some (or all) of the above mentioned effects are explained after all by some different 
physical (and geometrical) reasons, not related to the cosmic rotation. However, "whether 
our universe is rotating or not, it is of fundamental interest to understand the interrelation 
between rotation and other aspects of cosmological models as well as to understand the 
observational significance of an overall rotation" | |369 . 



We believe that the cosmic rotation is an important physical effect which should find its 
final place in cosmology. In this paper we outlined one of the possible theoretical frameworks 
which can underlie our understanding of this phenomenon. 
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